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Abstract 

We present a study of the QCD interactions which do not conserve the chirality of quarks. These 
non-perturbative forces are responsible for the violation of the Ua{^) charge conservation and for 
the breaking of chiral symmetry. From a systematic analysis we argue that the leading sources 
chirality flips are the interactions mediated by topological vacuum field fluctuations. We study in 
detail the contribution of instantons and derive a simple model-independent semi-classical predic¬ 
tion. This result can be used to check on the lattice if instantons are the dynamical mechanism 
responsible both for chiral symmetry breaking and for the anomalous violation of the Fa( 1) charge 
conservation. 
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I. INTRODUCTION AND MOTIVATION 


Understanding the structure of the quark-quark interaction at all scales is a fundamental 
task of nuclear and high-energy physics, which requires the solution of the non-perturbative 
sector of QCD. Although this goal is not yet been completely achieved, a good deal of 
information has been gathered in the past decades. In particular, we known that the non- 
perturbative dynamics mixes quark modes of different chirality. This property of the strong 
interaction is signaled by the anomalous violation of the Ua( 1 ) charge conservation and by 
the breaking of chiral symmetry. These two phenomena have important consequences on the 
physics of the light hadrons. The anomalous U^(l) Breaking (UIB) allows to explain the 
absence of a ninth Goldstone boson. The Spontaneous Chiral Symmetry Breaking (SCSB) 
shapes the structure of the spectrum of the lightest hadrons and their interaction, at low 
momenta. Clearly, identifying the dynamical mechanism responsible for the SCSB and the 
UIB is a fundamental step toward our comprehension of the quark-quark interaction. 

From the dynamical point of view, the structure of the spectrum of the lightest mesons 
indicates that quark-quark interaction is particularly attractive in the flavor non-singlet 0 “ 
channel and less attractive (or possibly, at some distance, even repulsive) in the singlet 0 “ 


channel [t]'). Very useful insight in the p 
the Nambu Jona-Lasinio (NJL) model 


Iiysics of chiral symmetry breaking has come from 
. In this approach, one postulates an effective 
quark-quark interaction, which breaks spontaneously the Uv{S) x 1/^(3) symmetry down to 
Uy(3). Such a symmetry breaking pattern would obviously generate nine Coldstone bosons, 
being also the U^(l) symmetry broken spontaneously. Hence, in order to reproduce the 
observed 77 — 77 ' splitting, an additional term is introduced, which breaks U^(l) explicitly 
and simulates the contribution from the’t Hooft determinant interaction. Clearly, in such 
a framework, there are two independent sources of violation of the Ua(1) charge, because 
the dynamical origin of the spontaneous Ua{^) breaking and of the explicit (i.e. anomalous) 
UIB are distinct. In fact, by suppressing the instanton contribution in the Lagrangian, the 
77 — 77 ' splitting disappears, yet chiral symmetry remains broken, hence the pions remain 
light. 

Conversely, one can imagine the opposite scenario in which both chiral symmetry breaking 
and the rj — p' splitting are consequences of the same dynamical mechanism. In other words, 
one can ask whether the same gauge configurations which generate the quark condensate 
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also break the Ua{^) symmetry through the axial anomaly. In such a scenario, the Ua{^) 
charge conservation is violated only through the anomaly. Moreover, there cannot be an 
7] — T]' splitting without the SCSB, and vice-versa. 

Instantons provide an example of gauge conhgurations which solve the U{1) problem 


a„d. at the sa^e thne. bteah ehita. syntntetty Q (fot a teview see also 0), The staad- 
ing question is whether these semi-classical helds are the dominant conhgurations in both 
phenomena, or if their contribution is only sub-leading. Unfortunately, since a systematic 
semi-classical approach to QCD is not possible, one has to rely on phenomenological models. 


B) 


such as the Instanton Liquid Model (ILM) l2| . This model was shown to quantitatively 
reproduce the breaking of chiral symmetry in the vacuum, and its restoration at high tem- 
)eratures. Moreover, it also provides a good description of the spectrum light hadrons 


IflL llll. I 12 L Il3j | and their form factors llM lla ll7| . 

On the other hand, a traditional argument against the hyp othesis that instantons drive 
both the UIB and the SCSB is based on the large Nc limit In this limit, the topological 
susceptibility disappears, instantons are suppressed by the exponent of the action, yet chiral 
symmetry remains broken. One is then lead to suppose that, if the large Nc world is at least 
qualitatively similar to the real one, instantons might not be the leading mechanism for the 
SCSB and that the rj — 7]' splitting and the SCSB have a different dynamical origin. 

In a recent work, Schafer studied in detail the instanton content of QCD with many 


Q. 


colors IlSll. He observed that the large entropy of instantons in SU{Nc) can overcome the 
exponential suppression due to the action, in such a way that the instanton density can 
remain hnite, at large Nc- From an analysis based on numerical simulations and mean-held 
estimates, he found that in the ILM the quark condensate is of 0{Nc), while the t]' mass is 
of 0{1/Nc). From this facts, he concluded that the ILM is not necessarily in conhict with 
the large Nc analysis. Although a number _of lattice studies seem to conhrm a picture in 
which instantons drive the SCSB (see e.g. 


on these issues has not yet been reached 


and references therein), a general consensus 
and further studies are therefore needed. 

In this work, we set-up a framework to study the time-evolution of the chirality of quarks, 
and we apply it to investigate the dynamics underlying the UIB and the SCSB, in QCD. 
In the next section, we shall dehne two functions of the Euclidean time, R^{t) and 
which measure the rate of chirality-flips, in a quark-antiquark system with zero total angular 
momentum, and total isospin 0 and 1, respectively. We will show that the typical time 
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between two chirality-flipping interactions scales with the inverse of the rj' mass, in the 
J = 0 channel, and with the inverse of the mass of the 5 meson (which is the axial partner 
of the pion) in the 1=1 channel. Moreover, by performing a spectral analysis of the 
our probability amplitude ratios, we shall show that chirality flips are dominated by the 
anomalous breaking of Ua{1), while the contribution of the SBCS is only sub-leading. This 
result is model independent, because it is based only on the numerical values of the masses 
of the lowest lying scalar and pseudo-scalar mesons, which are known experimentally. It 
implies that, in the \qq) system under consideration, the leading chirality mixing interaction 
is mediated by topologically charged gauge fields. 

Then, we will use the Operator Product Expansion (OPE) to perform a systematic study 
of the dynamical origin of the mixing of chirality in QCD. We shall see that, if vacuum 
held huctuations were small (vacuum dominance approximation), then the chirality hips 
would be dominated by the SCSB, and the contribution from the anomalous UIB would be 
completely negligible. However, this prediction is not supported by phenomenology which, 
as we have mentioned above, seems to favor the opposite scenario in which the anomalous 
UIB dominates. Therefore, one is lead to conclude that vacuum huctuations play a major 
role in the mixing of chirality. Collecting all these observations, we shall argue that the 
leading dynamical source of helicity hips is represented by the interaction mediated by some 
large topological vacuum huctuations, which badly violate the factorization assumption. 

In section uni we will address the question of whether such dynamics is driven by in- 
stantons. We shall hrst compute analytically the leading instanton contribution to the 
helicity-hip probability amplitudes, in the ILM. Then, we shall derive a simple semi-classical 
prediction, which does not depend on the phenomenological parameters of the ILM. Such a 
relationship is a model-independent signature of the instanton-induced interaction and can 
be used to check in an unambiguous way if instantons are responsible both for the SCSB 
and for the UIB. The feasibility of checking such a signature using lattice simulations will 
be discussed. 

Interestingly, we find that the single-instanton effects generate the same rate of chirality 
flips, in the singlet and non-singlet channel. The equality R^{t) = (which holds 

for small values of r, for which many-instanton effects can be neglected) is non-trivial and 
represents another model independent signature of instanton-induced forces. 

Results and conclusions of this work are summarized in section HU 
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II. MIXING OF QUARK CHIRALITY IN QCD 


For sake of simplicity, let us consider QCD with two flavors. Moreover, since we are 
interested in the helicity flips generated by the quark-quark interaction, it is convenient to 
work in the chiral limit. This way, the purely kinematical chirality flips induced by the 
quark mass are suppressed. 

We begin by dehning the following combination of gauge invariant correlation functions: 

:= {0\T[uiT)PLd{T)d{0)Pnui0)\0) + iPL^Pn) 

+ {0\T[d{T)PLu{T)u{0)PRd{0)\0) + {PL^PR) (1) 

lmT[Hr)PLuiT)di0)PRdi0)\0) + iPn^PL) 

+ (0|T[d(r)Pid(r) d{0)PRd{0)\0). + {PR^PL) 

+ (0| T[d(r) Pl d{T) h(0) Pru{0) |0) + {Pr ^ Pr) 

+ (0| T[u{t) Pl u{t) h(0) P,jn(0) |0). + {Pr ^ Pr)] (2) 


where. 


Pr- = 


1+75 

2 


Pl-.= 


1-75 

2 


(3) 


denote the probability amplitude for a flavor singlet (non-singlet) 
\q q) state to be found, after a time interval r, in a state in which the chirality of the 
quark and antiquark is flipped^. In QCD, even in the case in which = rrid = 0, quark 
states with different chirality can mix under time evolution. Hence, we expect such matrix 
elements to be in non-vanishing, in general. 

Similarly, we dehne: 

ATn-fUp{r) := (0|T[n(r)Pid(r)J(0)Pin(0)|0) + (Pi^P«) 

+ (0|T[d(r)Pin(r)h(0)Pid(0)|0) + (Pi^P^) (4) 

^Ln-fiipi^) ■= ^[(0| T[h(r) Pl M(r) (i(0) Pl d(0) |0) + (Pl ^ P/j) 

+ (0| T[d(r) PLd(r) J(0) PLd(O) |0). + {Pr ^ Pr) 

+ (0| T[ d(r) Pl d(r) h(0) Pr n(0) |0) + {Pr ^ Pr) 

+ (0|T[n(r)PLn(r)h(0)PLn(0)|0). + (PL^Pi?)], (5) 


^ In all formulas, analytic continuation to Euclidean time is assumed. 
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These functions measure the probability amplitude for the quark and antiquark not to have 
exchanged their chirality, after a time interval r. 

It is convenient to rewrite such matrix elements in terms of scalar and pseudo-scalar 


mesonic two-point correlation functions: 

= n.(r)-n,(r) (6) 

AZ-fUpir) = n,(r) + n.(r), (7) 

= n<,(r)-n^,(r) (8) 

^Ln-fiipi^) = n<,(r)+ n^,(r), (9) 

where the correlation functions are dehned as^: 

n,(T) = {o|j,(T)jt(o)|o) (10) 

n,-(T) = {0|J,-(r)4(0)|0) (11) 

nj(T) = {o|Ji(T)j](o)|o> (12) 

n,(T) = {o|j,(T)jt(o)|o). (13) 


The interpolating operators, exciting states with given (J, J^) quantum numbers, are de¬ 
hned as: 


Jnir) 

:= M(r)i75d(r) 

(14) 

Jsir) 

:= u{t) dij) 

(16) 

Jn' (t) 

■= (M(r)i75n(r)-hd(r)i75d(r)) 

(16) 

Ja{T) 

:= ^ {u{t) u{t) + d{T) d{T)) 

(17) 


At this point it is worth observing that the operators Jjr, Js, Jp'-, and can be transformed 
into each other by means of appropriate Ua{^) transformation (see Fig. Q). 

With these amplitudes, we can now construct two ratios which represent the probability 
amplitude to hnd the quarks in the hipped chirality state, relative to the amplitude for 

^ In this expression, have denoted with rj' an iso-singlet pseudo-scalar state, i.e. the SU{Nf = 2) corre¬ 
spondent of the \r]o) state. Moreover we can, without loss of generality, choose r > 0 and disregard the 
Euclidean time-ordering. 
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FIG. 1: Ua{^) transformations relating the Jjr, Js-, Jr)'^ and J„ operators. 


remaining in the same state. In the singlet and non-singlet channels, we have: 



(18) 


(19) 


In the chiral limit, these ratios carry information about the rate of chirality mixing interac¬ 


tions, in QCD. 

Let us study m and (P|) in the limit of large Euclidean times. In this case, one expects 
quarks to interact several times. Intuitively, there should exist a characteristic time scale 
r, determining how often chirality mixing interactions occur in the QCD vacuum. After 
having exchanged their chirality many times, quarks essentially “loose the memory” of what 
was their initial state, therefore they should be found in either chirality conhguration, with 
equal probability. That is to say: 




( 20 ) 


1 . 


In order to derive this result rigorously and to determine f, we use the fact that, at large 
enough r, the spectral representation of the correlation functions (irniiT!?!) is dominated by 
the contribution of the lowest lying states with the appropriate quantum numbers 0: 



( 21 ) 


( 22 ) 


(23) 


( 24 ) 
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where Xn,s,a'q' are the coupling constants of the densities (IT31IT7I) to the corresponding lowest 
lying states. 

Hence, in the large r limit, we can rewrite the ratio in the following way: 

1 - 2 5(r) + o(e), (25) 


where 




(r-oo) / 

n.(r) \xj 


(26) 


These relationships show that the information about the initial chirality of the quarks is 
exponentially destroyed, as the time increases. From the exponent in (|2ni) can immediately 
read-off the corresponding characteristic time scale: 


= --- ~ 0.2 fm. (27) 

{ms - m^) 

We shall now prove that the ratio R^^{t) represents the amplitude for the chirality-flips 
induced by the non-conservation of the Ua{1) charge. If the Ua{1) symmetry is not spoiled, 
then we have = ms and = A^. As a consequence ^{t) ^ 1 and the chirality-flip 
amplitude ratio (P|) tends to zero, at inhnity. This condition is sufficient to imply that 
R^^{t) is identically zero. In fact, let us suppose that this function is not identically zero, 
yet vanishes at inhnity. This corresponds to the unphysical scenario in which quarks can hip 
their chirality only in some hnite time interval, but eventually have to “go back” to their 
initial chirality state. 

The fact that n 7 r(r) ^ n 5 (r) and therefore R^^{t) ^ 0 only implies that the axial charge 
is not conserved, regardless if this is predominantly because of a spontaneous or anomalous 
breaking of the 77^(1) symmetry^. The combination of correlators that relates directly to 
the anomaly is the diherence between the tt and two-point functions, n 7 r(r) — n,j/(r). It 
is therefore convenient to add and subtract the n,^/(r) correlator, from the numerator in 
i?^5(r): 



n^(r) - n,^'(r) lin'ij) - n^(r) 
n^(r) + n5(r) n^(r) + n5(r) 


=: Rllm.ir) + RsuA^) 


NS 


(28) 


^ Indeed, in the NJL model, in which the Fa( 1) is broken spontaneously, one has > n 5 (T). We thank 

T.Schafer for pointing this out. 


8 









We are now in condition to disentangle the violation of the axial charge conservation induced 
by the anomalous and by the spontaneous breaking of Ua{X)- Indeed, if there was no 
anomaly, and the Ua{^) symmetry was just spontaneously broken, then n^(r) = n^/(r) 
and Ranom.{'^) would vanish. In this case, the contribution to would come entirely 

from Rsua^^)- Conversely, if the symmetry was not broken, then we would have 

(r) (see Fig. Q), therefore Rsua{^) — 0- In the real world we observe a relatively 
large tt—t]' splitting. As a consequence, R^^{t) turns out to be dominated by the anomalous 
term, Ranom.i'^)- To see this, we use again the spectral representation, in the large Euclidean 
time limit; 

Ranon^ir) = 1 - C^r^, g-K—C - _ " + ..., (29) 

e-K—-) - (1 _ ^ g-K'—O (39) 


where 



(31) 


and ellipses denote higher order terms in -r yy^g ggg that, in this limit, 

Ranom.i.'^) of Order 1, while Rsua^'^) I® exponentially suppressed by the rj' — n mass differ¬ 
ence. We conclude that the violation of the axial charge conservation, parametrized by the 
function R^^{t), is mainly due to the anomalous breaking of the Ua{^) symmetry. This 
conclusion has quite important implications on quarks and gluons dynamics. It suggests 
that, in this channel, the leading chirality mixing interaction is mediated by topologically 
charged gauge helds. 

For sake of completeness, let us now consider the helicity-flip amplitude ratio in the 
singlet channel; 


R^t) 


1 — ( t —^00) 

1 + e(r) ^ 


l-2e(r) + o(e2). 


(32) 


where. 


e r ; = 


(t— ^ri' y/RRf 


g-m,,/ T 


(33) 


n^(r) v/2 8(0|gg|0)2 (7rr)3/2- 

In this channel, we could not hnd a simple way to single-out the effects due to the anoma¬ 
lous UIB. This is essentially because the scalar density operator (0 has a non-vanishing 
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vacuum expectation value (Eq. (j2Sl))- As a consequence, R^{t) receives an explicit contri¬ 
bution from the quark condensate^ (Eq. (Idd|l i. Notice that the anomalous UIB contribution 
induce an exponential mixing of chirality, while the SCSB participates only through an ad¬ 
ditional pre-exponent. The characteristic time, which determines the exponential mixing of 
chirality is the inverse of the r]' mass: 


Notice that -C which suggests that the non-perturbative physics associated 

with the UIB and SCSB is characterized by an additional scale, signihcantly larger than 
Rqcd- 

The phenomenological analysis performed so far has indicated that the helicity flips in 
the non-singlet channel are dominated by the dynamics responsible for the anomalous UIB. 
Now we want investigate what can be learnt from this result about the non-perturbative 
quark-quark dynamics. At this purpose, we need to consider in the opposite limit 

of small Euclidean times, where quark and gluons are the relevant degrees of freedom. After 


performing Wick contractions, the correlation functions (fTUHTl read: 

n^(r) = (Tr [S'(r, 0 ) 75 S'( 0 ,r) 75 ]), (35) 

n 5 (r) = -(Tr [^(r,0)^(0,r)]), (36) 

n^/(r) = [(Tr [^(r, 0 ) 75 ^( 0 ,r) 75 ]) -2 (Tr [^(r, r) 75 ] Tr[^(0, 0) 75 ])], (37) 

n.(r) = [2 ( Tr [S{t, r)] Tr[^(0, 0)]) - (Tr [S{t, 0) ^(0, r)])], (38) 


where S'(r, 0) is the quark propagator and (■) denotes the average over all gauge conhgu- 
rations. At extremely small distances, asymptotic freedom demands that the correlators 
(unmnD approach those evaluated in the free theory. One has: 

n,(r),n,,(T),ni(T),n,(r) n„(r) := ( 39 ) 

71 * 

Consequently, we have R^{t), R^^{t) 0 . This is an expected result: in a free theory 

of massless quarks, chirality is a conserved quantum number, and the rate of quark helicity 
flips is zero. 


^ It is tempting to simply replace Ho-lr) with a modified scalar two-point functions, in which the vacuum 
contribution has been subtracted out, n(^(T) := ncr(T) —2(0|g g|0)^. However, the corresponding amplitude 
ratio, R'^{t) := (n(, — n,;/)/(n(,. -|- H,,/) would not have a simple probabilistic interpretation in terms of 
chirality flips. 
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As r gradually increases from zero, perturbative corrections to the correlators (ITUIIT^ 
become more and more important. However, their contribution to m and (USD has to 
vanish exactly, because the perturbative quark-gluon vertex does not mix quark modes with 
opposite chirality. Hence, any deviation from zero of the ratios (UHl) and (uni) is a signature 
of non-perturbative physics. 

Non-perturbative QCD corrections, in the limit r —0, can be evaluated in a systematic 
way, using the Operator Product Expansion (OPE). It known that, in the chiral limit, the 
quark and gluon condensate singular terms in OPE cannot distinguish between scalar and 


22 L Hence, including the contribution up to dimension 


pseudo-scalar mesonic correlators 
4 operators one hnds = 0. 

In the’t Hooft picture, in which the anomaly is realized by instantons, the fact that such 
singular terms cannot reproduce the chirality flipping due to the non-conservation of the axial 
charge is an expected result. In fact, the leading instanton contributions to the correlators 
(Unmi tend to a constant, in the r —0 limit. Therefore, one expects some deviation from 
zero in to appear, starting from the regular terms, i.e. from dimension 6 operators. 

Indeed, we hnd^: 

. -2 TT^ 


R^^ (t) ~ 

anom. \ ) 


Nr. 


-{Olui'y^udi'f^dlO) r° -|- .. 


jNS 


IT 


Rou Ir) ~ 

SUa^ ! JY 


(Ojuiys M uiys m|0) — {0\uuuu\0) 


(Oluiys M diys (i|0) -|- {0\uudd\0) 




(40) 


(41) 


Let us now discuss the implications of this result. Quite interestingly, we have found that 
the quantity: 


A : = lim 

r^O 


R^SAir) 


R^^ (t) 

^anom. \ ' ) 

(Oluzys u ufys m| 0) -|- (OlhfysudfysdlO) -|- {Q\uu dd\Q) — {Q\uuuu\Q) 


(42) 


2 (OluiysMdiysdlO) 

parametrizes the relative contribution to quark helicity flips of the spontaneous symmetry 
breaking with respect to the anomalous symmetry breaking. In particular, if A ~ 0, then 
the violation of the Ua{R) charge comes almost entirely from the anomaly. 


® In this expression and in (03 we have used isospin symmetry and we have neglected logarithmic pertur¬ 
bative corrections. Notice that the term proportional to the quark condensate and that proportional to 
the mixed quark condensate g{Q\qG^v<y^ivq.Wl are suppressed, in the chiral limit. 
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It is instructive to study the behavior of the chirality flip ratio in the large Nc 

limit. In this limit, the spectral representation of the condensates is dominated by vacuum 
insertions and one can use the factorization approximation: 

(0|g T,q q q\0) ~ ^(Tr[ri] TiiT^] - Tr[ri hs]) {0\qq\0y, Af:=4Nf N,, (43) 

which gives (Oluij^udij^dlO) ~ 0, hence A —> cx). This result implies that, in the asymp¬ 
totically large world, the chirality flips are dominated by the spontaneous chiral symmetry 
breaking, while the anomalous contribution disappears, in agreement with common wisdom. 

In the real world, the accuracy of the factorization approximation is questionable. In¬ 
deed, the analysis performed at large distances indicates that the spin flips are actually 
dominated by the anomaly and not by the SCSB (Eqs. and (jSni))- In general, the 

vacuum dominance assumption is expected to fail in in the presence of large vacuum held 
huctuations. Hence, our analysis suggests that in QCD the helicity hips are mainly induced 
by some large vacuum gauge held huctuations. As we mentioned above, such huctuations 
have to be topologically charged, as they induce the UIB though the anomaly. The most 
natural candidates are instantons. Their contribution will be analyzed in the next section. 


III. INSTANTON-INDUCED CHIRALITY FLIPS 


In the previous section we have argued that the helicity hips are dominated by some large 
held huctuations which have connection with the axial anomaly. Instantons are example of 
gauge helds which, at the same time, provide a realization of the anomaly and break chiral 
symmetry. It is therefore interesting to estimate their contribution to the ratio (pll. At 
short distances, the leading instanton contribution to the relevant correlation functions (HU 
can be computed analytically, using the Single Instanton Approximation (SIA). This is 
an ehective theory of the instanton vacuum (for a detailed description of this approach, see 
j^), in which the contribution of the closest instanton is taken into account explicitly, while 
all other (infra-red) multi-instanton degrees of freedom are integrated out and replaced by 
a single effective parameter m*, dehned as: 

1/2 

n \ 1 


m = 


dpd{p) 


Stt^ P^ 


TrE/,jV’o/(a:) 


(44) 
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where d{p) is the instanton size distribution, n is the instanton density, 'ipopx) is the zero 
mode wave function in the held of the J-th pseudo-particle and T/j is the overlap matrix, 
given by: 




(45) 


The effective parameter m*, computed numerically in the ILM from (m, is of the order of 
70 MeV. The main advantage of the SIA is that it allows to obtain the leading instanton 
contributions to Green’s functions by means of simple analytical calculations. The range of 
applicability of this approach was studied in detail 17|,|2^. It was found to be very accurate 
for correlation functions smaller than the typical distance between two neighbor instantons. 

The evaluation of the correlation functions leads to 


nf= n„(T) + no,(r) 

nf= n„(T) - 


(46) 

(47) 

(48) 


where no('r) is the usual free correlator dehned in ()d9|l and no 7 (r) is the instanton contri¬ 
bution, given by: 


Tloi{x) = 


An 


TT^ 


d!^z / dp d{p) 


P 


[z^ + p^f {{z - xf + P^Y^ 
The corresponding result for the chirality-hip amplitude ratio is: 

dns , _ no;(T) 

«s//i(T) - jj , , , 


(49) 


(50) 


In passing, we observe that the leading instanton contribution to R^{t) is the same as 
that to 


R 


s 

SIA 


{t) = R 


NS 

SIA 


[T . 


(51) 


This equality, which is valid only in the region in which the one-instanton ehects are dom¬ 
inant, is non-trivial. In fact, we recall that at large distances these quantities have very 


® For sake of simplicity, in this calculation we have replaced the non-zero mode part of the quark propagator 
with the free propagator (zero-mode approximation). It is possible to show that this approximation is 
very accurate for the particular correlation functions we are considering |2!tI| . The main results of this 
section (Eq. and Eg. (I54II below) do not change, when one correctly accounts for the non zero-mode 
contributions. 
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different spectral representations. Hence, in general, we expect them to be different func¬ 
tions of r. 

Let us now establish the connection with the OPE analysis of the previous section. Ex¬ 
panding Eq. (m for short Euclidean times one Ends: 

= / dpd{p)^T^ + ... (52) 

where the ellipsis denote terms which are higher orders in r. 

In general, calculations of instantonic effects require the knowledge of two quantities 
which cannot be obtained in a systematic way in QCD: the instanton density h, and the 
instanton size distribution d{py. From a phenomenological estimate, Shuryak suggested to 
usej^ h ~ 1 fm“^, and d{p) = (5(p—l/3fm). In Fig. 0 the SIA prediction (IKn|l for is 

compared with its short-time expansion (El, and with the result of numerical simulations in 
the Random Instanton Liquid Model (RILM)®, in which many-instantons effects are taken 
explicitly into account. We observe that the very simple SIA analytic prediction agrees 
with the much more complicated RILM calculation up to quite large times r ~ 0.6 fm. 
Moreover, we observe that the short-time expression (El converges up to distances of the 
order of r < 0.25 fm. This provides an estimate of the radius of convergence of OPE 
expansions, when the direct-instanton term is included. 

The results discussed so far depend on the phenomenological parameters of the ILM. 
Now we make one more step and show that it is possible to circumvent all the model 
dependence®, and obtain the semi-classical prediction. At this purpose, we observe that the 
instanton contribution to the four-quark condensate reads: 

2 n f 1 

{yqqqq\Q)siA = ^ / dpd{p)—= y)\q-f^ qq-f^q\Qi )sia (53) 

5 m ^ J p^ 

We notice that leading instanton effects contribute equally to the pseudo-scalar and scalar 
four-condensates and therefore badly violate the factorization assumption. This implies 
that, in the instanton vacuum, A ~ 0. Therefore, the instanton-induced UIB is channeled 

^ We recall that the effective parameter m* can be obtained from d{p) and n through Eq. 14411 . 

® For a detailed description of this model and a comparison with other versions of the ILM, see |9|. 

® As another example of model-independent instanton prediction in QCD is discussed in [ij. In that paper, 
we computed the proton and pion dispersion curves and we derived, in a parameter-free way, the instanton 
contribution to the nucleon mass. 
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FIG. 2: The chirality flipping amplitude ratio, in the non-singlet channel, Points are 

numerical simulations in the RILM, the solid line denotes the SIA prediction (EDI), and the dashed 
line its short-time expression •The dot-dashed line represents the explicit breaking contribution 
(jnSJ, with m = 150 MeV. 

exclusively through the axial anomaly, in agreement with the results of the phenomenological 
analysis presented in the previous section. 

Combining dsni) with (ESI we can include all unknown quantities appearing in our calcu¬ 
lation in the expression for the quark condensate. This way, we obtain the direct-instanton 
contribution to OPE for m and (Hi: 

= y + (54) 

It is worth emphasizing that this semi-classical calculation requires no assumption on the 
instanton size and density, yet predicts a particular relationship between different vacuum 
expectation values. Hence, if the semi-classical contribution is large, so that instantons are 
the leading conhgurations responsible for the UIB and the SCSB, then this relationship 
must be at least approximatively satished, also when the quantum average is performed 
over all conhgurations (and not only over the semi-classical huctuations). 

Lattice QCD represents the natural framework in which performing such a test. On the 
practical level, however, up-to-date simulations still need to use quite large values of the 
current quark masses. Far from the chiral limit, the explicit breaking of chiral symmetry 
provides a competing source of chirality hips, which has to taken into account. The impor- 
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tance of this effect can be estimated by evaluating (HI in the free theory - using a massive 
quark propagator - and then comparing with the ILM prediction dSl. After an elementary 
calculation, we obtain: 

R^Ucu{r) = ^r^ + ... (55) 

where m is the current quark mass. Clearly, at short distances, the contribution to helicity 
flips coming from the explicit breaking of chiral symmetry dSSl) will necessarily dominate 
over the dynamical effects discussed above. 

The key question is if there exists a window of Euclidean times in which the instanton 
effects are dominant over the kinematic ones, and the single-instanton prediction (P|) is 
reliable. In Fig. m the contribution of explicit chiral symmetry breaking (HU) with a current 
mass m = 150 MeV (of the order of those used in a typical lattice calculation) is compared 
with the SIA and RILM estimates. We conclude that, in the region 0.1 fm < r < 0.25 fm, 
the instanton contribution should be dominant. Moreover, in this window, the relevant 
correlation functions (Hnma) are of order 1, so it should not be numerically very challenging 
to measure this signal on the lattice. 


IV. CONCLUSIONS 

In this paper we have presented a study of the mixing of chirality under time evolution, 
in QCD. In the chiral limit, quark helicity-flips have a purely dynamical origin, therefore 
they carry information about the quark-quark interaction. 

We have constructed two combinations of correlation functions, and which 

are related to the rate of chirality flips in a quark-antiquark pair, propagating in the QCD 
vacuum. Such functions receive no contribution from any number of perturbative gluon ex¬ 
changes, therefore they represent useful tools to investigate the non-perturbative dynamics. 

We have shown that the chirality flipping processes contributing to R^^{t) are induced 
by the same non-perturbative forces responsible for the non-conservation of the axial charge. 
In fact, the ratio R^^{t) would be identically zero, if the degeneracy between the pion and 
its Ua{1) partner was not lifted. Using the spectral decomposition at large Euclidean times, 
we have studied the relative contribution to the UIB coming from the axial anomaly and 
from the possible spontaneous breaking of ^ 4 ( 1 ). We found that the UIB comes almost 
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entirely from the axial anomaly. From a such a spectral analysis, we have also obtained 
an estimate of the characteristic time between two consecutive spin-flipping interactions 
in QCD and found f ~ 0.2 fm. This value suggests that the chirality mixing dynamics 
is characterized by a non-perturbative length scale which is signihcantly smaller than the 
typical conhnement scale, 1/Kqcd — 1 

We have calculated the first non-perturbative corrections to the amplitude ratio 
using OPE. We have found that these are associated to the regular terms, which are com¬ 
monly related to the so-called “direct instanton” contribution. This is consistent with the 
’t Hooft solution of the U{1) problem, in which the axial anomaly is realized by instantons. 
We observed that the fact that is dominated by the anomaly implies that vacuum 

held huctuations are large. In fact, in the vacuum dominance approximation, in which 
huctuations are neglected, the chirality hips are dominated by the SCSB, in disagreement 
with what we found phenomenologically from the spectral analysis. This suggests that the 
helicity hips are induced by topological gauge held huctuations. 

We calculated analytically the instanton contribution to R^{t) and R^^{t) in the ILM, 
using the SIA. Then, we derived a model-independent semi-classical relation, which connects 
scalar and pseudo-scalar correlators to the four-quark condensates. This prediction repre¬ 
sents a clean signature of instanton-induced dynamics. We suggest that it should be checked 
with an ab initio calculation, on the lattice. We discussed the feasibility of performing such 
a calculation and concluded that the artifacts related to large quark mass ehects should be 
negligible. 

An application of the same framework to the study the quenching ehects in lattice QCD 
is in preparation. 
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